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Abstract 

In the heavy quark limit, the heavy baryon A q (Q=b or c) can be regarded 
as composed of a heavy quark and a scalar light diquark which has good spin 
and flavor quantum numbers. Based on this picture we establish the Bethe- 
Salpeter (B-S) equation for Aq in the leading order of 1 /tuq expansion. With 
the kernel containing both the scalar confinement and one-gluon-exchange 
terms we solve the B-S equation numerically. The Isgur-Wise function for 
Ab —> A c is obtained numerically from our model. Comparison with other 
model calculations are also presented. It seems that the Isgur-Wise function 
for At, —> A c drops faster than that for B —> D. The differential and total 
decay widths for A& —> A c lP are given in the limit rrib,c oo. 


PACS Numbers: 11.10.St, 12.39.Hg, 14.40.Nd, 14.40.Lb 


I. Introduction 


The past few years have seen much progress in heavy flavor physics due to the 
discovery of the new flavor and spin symmetries SU (2) f x SU (2) s in the heavy quark 
limit and the establishment of the heavy quark effective theory (HQET) [[T]. In the 
framework of the HQET the processes involving heavy quarks can be simplified 
since the HQET may reduce the number of form factors. For instance in the leading 
order of 1/rriQ expansion only one form factor (the Isgur-Wise function) remains for 
A b —> A c transition. However, to get the complete physics one still has to face the 
basic problem, the nonperturbative effects of QCD. As the formally exact equation 
to describe the relativistic bound state, the B-S equation @0 can be simplified 
to a great extent in the heavy quark limit and has been applied to give many 
theoretical results concerning heavy mesons Q. At present, there are not many 
experimental results for heavy baryons. But in the future we may expect more and 
more data coming out in LEP and the forthcoming B-factory. Actually recently 
OPAL has measured some physical quantities for A;, such as the lifetime of A;, and 
the product branching ratio for Aft —> Al~uX [j5|. Hence the study of heavy baryons 
are of particular interestes in the near future. It is the motivation of the present 
paper to establish the B-S equation for the heavy baryon in the heavy quark limit 
and then solve it numerically by assuming some reasonable form of its kernel, and 
consequently give some phenomenological predictions. 

When the quark mass is very heavy comparing with the QCD scale A qcd, the 
light degrees of freedom in a heavy baryon A q (Q=b or c) becomes blind to the 
flavor and spin quantum numbers of the heavy quark because of the SU (2) f x SU (2) s 
symmetries. Therefore, the angular momentum and flavor quantum numbers of the 
light degrees of freedom (the light cliquark) become good quantum numbers which 
can be used to classify heavy baryons. For example, A q and Eq correspond to 
the zero and one angular momentum states of the light diquark respectively. Hence 
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it is natural to regard the heavy baryon to be composed of a heavy quark and a 
light diquark. In the present paper we concentrate on Aq which has a scalar light 
diquark with [ud] flavor quantum number and zero spin and isospin. Other 
heavy baryons like can be dealt with in the similar way. 

Based on the above picture of the composition of the heavy baryon the three 
body system is simplified to two body system. We will establish the B-S equation 
for the heavy baryon in this picture. Then we will solve this equation numerically by 
assuming that the kernel contains the scalar confinement and one-gluon-exchange 
terms. Then we will apply the B-S equation to calculate the Isgur-Wise function for 
A b —■> A c and the decay width for the semileptonic decay of A b to A c in the order 
m b , c -> cx). 

The remainder of this paper is organized as follows. In Sect. II we establish 
the B-S equation for the heavy quark and light scalar diquark system and discuss 
the form of its kernel. Consequently this equation can be solved numerically. Then 
in Sect. Ill we apply the solution of the B-S equation to calculate the Isgur-Wise 
function for A b —> A c and also the differential and total decay widths for A;, —> A c li>. 
Comparison with other model calculations for the Isgur-Wise function will also be 
presented. Finally, Sect. VI is served for summary and discussions. 

II. The B-S equation for A q 

As discussed in Introduction Aq is regarded as the bound state of a heavy quark 
and a light scalar diquark. Based on this picture we can define the B-S wave function 
of Aq by 

x(xi,x 2 ,P) =< 0 \Ti/)(xi)<p(x 2 )\Aq(P) >, (1) 

where ^(xi) and (p(x 2 ) are the wave fuctions of the heavy quark and the light scalar 
diquark repectively, P = 'm\ Q v is the momentum of Aq and v is its velocity. Let 
rriQ and mo be the masses of the heavy quark and the light diquark in the baryon, 
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Ai = ——, A 2 = — ,J1 r —1 V be the relative momentum of the two constituents. 

1 rriQ+rriD ’ mQ+mr> ’ r 

The B-S wave function in the momentum space is defined as 

x(xi, x 2 , p ) = e lPX J e tpx xp(p), ( 2 ) 


where X = Aiaq + X 2 x 2 is the coordinate of the center of mass and x = Xi — x 2 . 
The momentum of the heavy quark is pi = X\P + p and that of the diquark is 
p 2 = -A 2 P + p. 

The derivation of the B-S equation for the two fermion systems can be found in 
textbooks |3|. In the same way we can prove that for the fermion and scalar object 
system the form of the B-S equation is still valid. Xp('P) satisfies the following B-S 
equation 

Xp(p) = Sf(AiP + p) j ■0^G(P,p,q) X p(q)S D (-\ 2 P +p), (3) 


where G(P,p,q) is the kernel which is defined as the sum of the two particle irre¬ 
ducible diagrams. In the following we will use the variables pi = v ■ p — X 2 m \ Q , p t = 
p — (v ■ p)v. Then in the leading order of 1 /itiq expansion we have 


Sp{X\P + p) 
Sd(~X 2 P + p) 


i(l + /) 

2 {pi + E 0 + m D + ie) ’ 
i 

pf-Wt + Eo + ic’ 


where W p = \Jpt + m 2 D an d E 0 is the binding energy. The corrections to eq. 
from 0(1/mo). 

In general, Xp(p) can be expanded as the following form, 


( 4 ) 


are 


Xp(p) = ( A + B f + Cf + D/f)u Kq {v , s), (5) 


where u\ Q (v, s) is the spinor of A q with helicity s and A, B, C, D are Lorentze scalar 
functions. 
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From eqs. (|3|) and ([|) it can be seen that 


t’Xp(p) = Xp(p)- ( 6 ) 

Combining eqs. @ and (|j) we immediately see that 

Xp(p) = 0p(p)ua q (v, s), ( 7 ) 

where 4>p(p) is a scalar field. This form is consistent with our picture about the 
heavy baryon. 

We assume the kernel has the form 

-iG = I ® IV i + ® (p 2 + P^)^, ( 8 ) 


where the first term arises from scalar confinement and the second one is from one 
gluon exchange diagram. The vertex of gluon with two scalar diquarks is propor¬ 
tional to (p 2 + p' 2 )fi {p 2 and p’ 2 are the momenta of the two scalar diquarks on the 
vertex) and a form factor F((p 2 + P 2 ) 2 ) (see Fig. 1). Substituting eqs. (|]), (0) and 
(D into (|) we have 


<Mp) 


i 

C Pi + E 0 + m D + ie){jpi - W p 2 + ie) 


d 4 q 


G(P,p,q)<i>p(q). 


( 9 ) 


Now consider the vertex of two heavy quarks with gluon. The momenta of the 
two heavy quarks are pi = \irriA Q v + p and p\ = \imx Q v + q respectively, p and q 
are of the order A qcd- In the heavy quark limit the heavy quark is almost on-shell, 
therefore, pi = qi on this vertex. Hence we can make the convariant instantaneous 
approximation J4| in the kernel. Then eq. (|9]) becomes 


<Mp) 


i 

0 Pi + Eo + me + ie) (pf - !Hp 2 + ie) 


^ 4(^1 + 2 PiV 2 )(j)p{q), 


( 10 ) 


where V is defined as V = V\ Pl=qi 
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In general 4>p(p) can be the function of pi and p t . Defining 4>p(pt) = / ^<Ap(p) 


one gets immediately the B-S equation for (j)p(pt ) 

1 f d 3 q t 


<f>p{pt) = 


2 ( E 0 — W p + m D )W p J (27r) 3 


(Vi — 2 W p V 2 )Mq t ). 


( 11 ) 


If one knows the form for the kernel V\ and V 2 , then (j)p(pt) can be obtained. 
Consequently <f>p(p) can be solved since after integrating q t in eq. m we get the 
relation between dpijj) and <j>p(qt) 

i f d 3 q t 


<Mp) = 


C Pi + Eo + m D + ie)(pf - ID 2 + ie) J (2vr) 3 


(i i + 2 P iV 2 )Mq t )- (12) 


The kernel V\ and V 2 for the B-S equation in the meson case was given in [H] as 


14 


8lThi' 


[(Pt - qt ) 2 + p 


212 


(2^f5\p t - q t ) 


d 3 k 871 -k 1 

(27r) 3 (k 2 + /i 2 ) 2 ’ 


14 


167T 


2 meson 


a seff 


(13) 


(14) 


3 (pt - qt) 2 + /i 2 ’ 

where k’ and oi se ff are coupling parameters related to scalar confinement and one- 
gluon-exchange diagram respectively. The second term in eq. (|T3J) is the counter 
term which removes the infra-red divergence in the integral equation. From potential 
model k' is from 0.18 to 0.2. The parameter /j is introduced to avoid the infra-red 
divergence in numerical calculations. The limit /x —^ 0 is taken in the end. 

In the baryon case, since confinement is still due to scalar interaction the form 
of V\ need not be changed. Only the parameter k' has to be replaced by k which 
describes the confinement interaction between the quark and diquark. However, 
since the cliquark is not a point-like object, there should be a form factor F(Q 2 )(Q = 
P '2 ~P 2 ) 1° describe the vertex of the gluon and two scalar diquarks. In general, this 
vertex is (Fig. 1) 

igsk^(P2+P 2 YF(Q 2 ) 
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with A’s being the Gell-Mann color matrices and g s the strong coupling constant. 
The form of F(Q 2 ) was assumed as 0 


HQ 2 ) 


a seffQp 

Q 2 + QV 


(15) 


where Q q is a parameter which freezes F{Q 2 ) when Q 2 is very small. In the high 
energy region the form factor is proportional to 1 /Q 2 which is consistent with per¬ 
turbative QCD calculations [0. By analyzing the eletromagnetic form factor for 
the proton it was found that Q q = 3.2GeV 2 can lead to consistent results with the 
experimental data [|5J]. 

Based on the above analysis the kernel for the B-S equation in the baryon case 
is of the following form 


V\ 


V 2 


87 TK 

[{jpt - qt) 2 + p 2 }' 2 


(27T) 3 5 3 (p t - q t ) 


d 3 k 87 tk 

(2tt) 3 {k 2 + fi 2 ) 2 ’ 


167T a 1effQo 

3 [{pt - Qt) 2 + h 2 ] [{pt - qt) 2 + Qo)' 


(16) 


There are two parameters k and a se ff in the kernel. However, there should some 
relation between them since when we solve the B-S equation (|TT|) numerically the 
binding energy should satisfy the following relation 


m AQ = m Q + m D + E (h 


(17) 


where we have omitted corrections from 0[l/mo) since we are working in the heavy 
quark limit. From the B-S equation solutions in meson case it has been found that 
the values mb = 5.02GeV and m c = 1.58GeV give predictions which are in good 
agreement with experiments 0 • Hence in the baryon case we expect 

mo + E 0 = 0.62GeH. (18) 


On the other hand the dimesion of k is three and that of k' is two. This ex¬ 
tra dimesion in k should be caused by nonperturbative diagrams which include the 
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frozen form factor F{Q 2 ) at low momentum region. Since A qcd is the only param¬ 
eter which is related to confinement we expect that 

k ~ A qcdk'. (19) 


Therefore in our numerical calculations we let k vary in the region between 
0.02GeV 3 to O.lGeV 3 . The diquark mass rrio is chosen to vary from 650 MeV to 
800 MeV. In order to satisfy the relation (p~8D we obtain the parameters for different 
values of mo- The results for mo = 700MeV are shown in Table 1. 


Table 1. Values of k and a se ff for m n — 700MeV 


^(GeV 3 ) 

0.02 

0.04 

0.06 

0.08 

0.1 

a seff 

0.67 

0.70 

0.72 

0.74 

0.75 


When m D = 650MeV we find that a se // = 0.60 for k = 0.02GeV 3 and « se // = 
0.71 for k = O.lGeV 3 . When m D = 800MeV, a se // = 0.79 for k = 0.02GeV 3 and 
ctseff — 0.84 for k = O.lGeV 3 . Having these two parameters we obtain the numerical 
solution for the B-S wave function directly. 


III. The Isgur-Wise function for —>• A c 


In this section we will apply the B-S equation for heavy baryons to obtain the 
numerical result for the Isgur-Wise function for A b —> A c . The transition diagram is 
plotted in Fig.2. In the limit m& )C —> oo there is only one form factor, the Isgur-Wise 
function to describe the weak transition from Af, to A c , 

< A c (n')|c 7 M 6 |A 6 (n) >= (20) 


where u — v ■ v' is the velocity transfer and U\ b and U\ c are the Dirac spinors of A b 
and A c respectively. 
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On the other hand the transition matrix element of A& —> A c is related to the 


B-S wave functions of A & and A c by the following equation 

< A c (n')|c7 /t 6|A ft (n) >= J 1 Xp'(p')Xp{p)Sd 1 (p2 ), (21) 

where P ( P' ) is the momentum of A*, (A c ). Xp'{p') is the wave function of the final 
state A c (i/) which satishes the constraint 

Xp'(p)/' = Xp'(p)- (22) 


The scalar part of the final state B-S wave function obeys the same B-S equation as 
(|9|). Then from eqs. (|20j) and (^1[) one gets immediately 

= / -^i ( t>P'{p)(t>p{p)SD 1 {P2)- (23) 

Since in the weak transition the diquark acts as a spectator its mementum in 
the initial and final baryons should be the same, p 2 = p 2 - Then we can show that 

p' = p + mp>(v' — v), (24) 


where again we omitted the 0(l/mg) corrections. From eq. ( |24|) we can get the 
relations between p' h p\ and pi , p t straightforwardly 


p'i = piuj - pt'/uj 2 - 1 cos9, 

p't = p^ + p^ (to 2 — 1 )cos 2 0 + p 2 (to 2 — 1) — 2pip t uVuj 2 — 1 cosd, (25) 

where 9 is the angle between p t and v[. 

Substituting the relation between 4>p(p ) and <t>p(pt) (Q) into eq. (|23|) and after 
integrating the pi component by selecting the proper contour we have 
d 3 p t 1 


pr,A = f d Pt _ _ _ 

“ J (2vr) 3 2W P (E 0 +m D - W P )(E 0 + m D - toW p - pW^ 2 - 1 cos6) 

/ dO f _ ___ _ 

j^^[ v i(Pt ~ r t ) - 2 (uWp + ptVu 2 - lcos9)V 2 {p' t - r t )\ \ k ,=-w P 4>p>{r t ) 
:[Vi(pt ~ k) ~ 2uW p V 2 ( Pt - l t )]Ml t ), (26) 


d 3 r f 
( 2 ? 
d 3 / t 

W) 













where the relation ( |25| ) has been substituted into the above equation. 


In eq. (p6|) all the three-dimensional integral can be simplified to one-dimensional 
integral. Furthermore, the normalization constant of the B-S wave function should 
be chosen such that the Isgur-Wise function is equal to one at zero momentum 
transfer. In Fig. 3 we plot numerical result for the Isgur-Wise function with the 
parameter mn = 700 MeV. For other values of rnp in the range of 650MeV to 
800MeV the shape of the Isgur-Wise function does not change a lot. 

The slope of the Isgur-Wise function at u — 1 which is of particular interests is 
defined as 


d£H 

do; 


0=1— —p~- 


(27) 


In Table 2 we list the values of p 2 for different rri£, and k. The values of p 2 
outside (inside) the braket correspond to n = 0.02GeV 3 (O.lGeV 3 ) respectively. 


Table 2. Values of p 2 


m D (GeV ) 

0.65 

0.7 

0.8 

P 2 

1.4 (2.4) 

1.4 (2.4) 

1.6 (2.4) 


It can be seen from Table 2 that the slope is insensitive to the values of mu. 
The Isgur-Wise function has also been calculated in other models. In Ref. || 
Guo and Kroll use the Drcll-Yan type overlap integrals for the model hadronic wave 
functions of and A c to obtain the following form of the Isgur-Wise function 


£M = 


u + 1 


exp 


-2 e z b 


2r 2 L 


uj — 1\ K 6 (2eb/yju + 1) 


u + 1 


K e (V2sb) 


(28) 


where e is the light scalar diquark mass and b is related to the average transverse 
momentum of the constituents in the heavy baryon. A") is defined as 


K t (x) — f dze z2 (z + x) 1 . 

J — X 
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The Isgur-Wise function obtained by Jenkins, Manohar and Wise [Of from the soliton 
model has the following form 

f(u;) = 0.99exp[-1.3(u - 1)]. (29) 


The MIT bag model calculation by Sadzikowski and Zalewski [|lT|) gives the 
following result 


= ( 


U + 1 


3.5+1.2/w 


(30) 


The values of p 2 from the above models are listed in the following 

2.9(3.7)[8] for b = 1.18(1.77 )GeV~\ 

P‘ = < 1.30. 


(31) 
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We can see that different models may give results for the Isgur-Wise function 
for A;, —> A c with some difference, p 2 > 1.3 in all these model calculations. The p 
parameter of the Isgur-Wise function for B —> D has also been calculated. In the 
B-S equation approach || p 2 varies from 1.0 to 1.1. Other model calculations for p in 


B —> D transition gives the following values: p 2 = 1.19±0.25[jnj]; p = 1.13±0.11||12 


and p = 1.20±0.17|T3|- Because of the large uncertainties we can not draw a definite 
conclusion by comparing the slopes of the Isgur-Wise fuction at the zero recoil point 
in the meson and baryon cases. However, it seems that the Isgur-Wise function drops 
faster in the baryon case than in the meson case. At least in the B-S approach the 
p parameter is bigger in the baryon case than in the meson case. 

The decay width for the semileptonic transition A & —> A c lu can also be predicted 
by using the Isgur-Wise function obtained. The differential decay width has the 
following form § [b| 


ur o 

— = -m\m Al A( 2 (u))Vui i ~ l[3u(i; + ?f‘) - 2 - 4u 2 ], 


(32) 


where p = m\ b /mf Ac and A = ^^ 3 1 V c b\ 2 B(A c —> ab ). \V c b\ is the Kobayashi-Maskawa 


matrix element. B(A C —> ab) is the branching ratio for the decay A c —> a(| + ) + 6(0 ) 
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through which A c is detected since the structure for such decay is already well known. 
The plot for is shown in Fig. 4 for mo = 700MeV. Again for other values of 

rrifj the results change only a little. 

After integrating u in eq. (fi2| ) we have the total decay width for A^ —»■ A c lP. 
For mo = 700MeV 


r(A b -> A c lv) 


4.2H(A C —> ab) x 10 10 s 1 when k = 0.02GeH 3 , 
5.7B(A C —> ab) x 10 10 s -1 when k = 0.1 GeV 3 . 


(33) 


When uid = 650MeV, r(A& —> A Ju) is 4.3(5.9)5(A C —■> ab) x 10 10 s 1 for k = 


0.02(0.1) and when m D = 800MeV, T(A b —> AJu) is 4.1(5.0)5(A C —> ab) x 10 10 s 1 


for k = 0 .02(0.1). 


IV. Summary and discussions 


In this paper we established the B-S equation for the heavy baryon which is 
considered as composed of a heavy quark and a light scalar diquark. By assuming 
the kernel containing scalar confinement and one-gluon-exchange terms we solve the 
B-S equation numerically. Furthermore, we applied the obtained results to calculate 
the Isgur-Wise function for A& —> A c . It is found that in the B-S approach, the 
Isgur-Wise function for A b —> A c decreeses faster than for B —>• D. Comparison 
with other model calculations are also presented. The slope parameter p 2 at the 
zero recoil point for At, —> A c is larger than 1.3 in all the present model calculations. 
In the limit m b)C —>• oo we also calculated the differential and total decay widths for 
A b —> A c lP. 

Although the B-S equation is formally the exact equation to describe the bound 
state, there is much difficulty in applying it to the real physical state. The most 
difficult point is the that we can not solve out the form of the kernel. Hence we have 
to use some phenomenological kernel which is reasonable because of the success 
of the potential model. This leads to some uncertainties. In our approach, the 
parameters k and ct S eff hi the kernel are not exactly determined. Therefore, we let 
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them to vary in some resonable range. The experimental data concerning A q in the 
future can help to fix the parameters in our model. 

We have worked in the heavy quark limit. The physical predictions such as 
T(A fe —> A c lv) will be subjective to 1/rriQ corrections. The study on the I/tttq 
corrections will appear in our later work. 
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Figure Captions 


Fig.l The diquark-gluon-diquark vertex. 


Fig.2 The weak transition diagram for A& —> A c . 

Fig.3 The Isgnr-Wise function for A& —> A c (rrio = 700MeV). For the solid line 
k = 0.02GeV 3 and for the dotted line k = O.lGeV 3 . 

Fig.4 The differential decay width for A;, —> A c lu (mo = 700MeV). For the solid 
line k = 0.02GeV 3 and for the dotted line k = O.lGeV 3 . 
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